function CrackIter
global glb
da=glb.crack.da;
N=glb.crack.N;

dN=100;%min([glb.geom.W/max(da)*0.5,1000]);

glb.crack.dN=dN;
glb.crack.N=N+dN;